High field domain wall propagation velocity in magnetic nanowires 
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A thory of field-induced domain wall (DW) propagation is developed. The theory not only explains 
why a DW in a defect-free nanowire must propagate at a finite velocity, but also provides a proper 
definition of DW propagation velocity. This definition, valid for an arbitrary DW structure, allows 
one to compute the instantaneous DW velocity in a meaningful way even when the DW is not 
moving as a rigid body. A new velocity-field formula beyond the Walker breakdown field, which is 
in excellent agreement with both experiments and numerical simulations, is derived. 



It is a textbook knowledge [l[ that a magnetic field can 
drive a magnetic domain wall (DW) to move. However, 
our understanding of the field-induced DW motion is far 
from complete although it has been intensively studied 
for more than fifty years and many interesting phenom- 
ena of magnetization dynamics have been found. Recent 
development in nanomagnetismQ demands a deep un- 
derstanding of DW motion in nanowires, especially how 
a field affects DW propagation velocity. DW dynamics 
is governed by the Landau-Lifshtiz-Gilbert (LLG) equa- 
tion that can only be solved analytically for some spe- 
cial problems 0, 0]. A number of theories have been 
widely accepted and written in books [H, such as kinetic 
potential approach that assumes zero-damping, Thielc 
dynamic force equilibrium formulation that is correct for 
rigid DW propagation, Schrycr and Walker analytical so- 
lution that is valid only for ID and exact only for field 
smaller than a so-called Walker breakdown field -ffivQj 
Slonczewski formulation that simplifies a DW by its cen- 
ter and the cant angle of DW plane. None of these or- 
thodox theories works beyond Hw although they have 
greatly enriched our current understanding of DW dy- 
namics. For example, kinetic potential approach cannot 
be a correct description of DW propagation because it vi- 
olates the principle of "no damping, no propagation" that 
will be explained in this paper. Thiele approach is a good 
way to describe a rigid DW propagation for small field 
H < Hw, but its assumptions are not valid for H > Hw- 
Schryer and Walker's approach is for ID and H < Hw, 
and its predictions for H > Hw are incorrect. For in- 
stance, its prediction that the v — H line for H >> Hw 
passing through the origin differs from both experiments 
and micromagnetic simulations [f| 0JE M- Its general- 
ization predicts a saturated velocity[5| (bounded by the 
velocity at Hw) that does not agree either with exper- 
iments or with simulations [1, 0> la 0- Slonczewski for- 
mulation is a great simplification of LLG equation that 
not only replaces partial differential equations by ordi- 
nary differential ones, but also is based on Thiele rigid 
DW approximation although the Slonczewski equations 
have also been applied to the case of H > Hw where it is 
known that DW deformation cannot be neglected. The 
problems with both Thiele and Slonczewski formulations 
can also be seen from their v — H formulaH G3 that 



do not capture the trend for H > Hw- Even more sur- 
prising, none of the existing theories provides a proper 
definition of DW propagation velocity when a DW does 
not propagate like a rigid body. 

In this paper, we develop a general formulation of DW 
propagation for both H < Hw and H > Hw that does 
not have all the problems with the existing theories. The 
theory reveals the origin of DW propagation. Firstly, 
we show that no static tail-to-tail (TT) or head-to-head 
(HH) DW is allowed in a homogeneous nanowire in the 
presence of an external magnetic field. A moving DW 
must dissipate energy because of various damping mecha- 
nisms. The energy loss should be supplied by the Zccman 
energy released from the DW propagation. This energy 
consideration can clearly explain DW velocity oscillation 
for H > Hw ■ Secondly, the energy conservation provides 
a proper definition for DW propagation velocity. This 
definition leads to a general relationship between DW 
propagation velocity and the DW structure. Finally, a 
new velocity-field formula beyond the Walker breakdown 
field is derived and is compared with both experiments 
and numerical simulations. 




FIG. 1: Schematic diagram of a HH DW of width A in a 
magnetic nanowire of cross-section A. The wire consists of 
three phases: two domains and one DW. The magnetization 
in domains I and II is along +z-direction (9 — 0) and -z- 
direction (6 — n), respectively. Ill is the DW region whose 
magnetization structure could be very complicate. H is an 
external field along +z-direction. 

In a magnetic material, magnetic domains are formed 
in order to minimize the stray field energy. A DW that 
separates two domains is defined by the balance between 
the exchange energy and the magnetic anisotropy energy. 
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To describe a HH DW in a magnetic nanowire, let us con- 
sider a wire with its easy-axis along the wire axis which 
is chosen as the z-axis as illustrated in Fig. 1. Since 
the magnitude of the magnetization M does not change 
according to the LLG equation^, the magnetic state of 
the wire can be conveniently described by the polar an- 
gle 9(x, t) (angle between M and the z-axis) and the az- 
imuthal angle <fi(x,t). The magnetization energy of the 
wire energy can be written in general as 



(1) 



E = J F{0,(t),Ve,V4>)d z x, 

F = f(6,<p) + - [(V0) 2 + sin 2 6>(V0) 2 ] - MH cos 9, 

where / is the energy density due to all kinds of magnetic 
anisotropics which has two equal minima at 9 — and 
7r (f(9 = 0, <f>) = f(9 = it, (f))). J describes the exchange 
energy. M is the magnitude of magnetization, and H is 
the external magnetic field along z-axis. In the absence 
of H, a static HH DW that separates 9 — and 9 = it 
domains can exist [3j in the wire. 

Non-existence of a static HH (TT) DW in a magnetic 
field-ln order to show that no intrinsic static HH DW 
is allowed in the presence of an external field (H ^ 0), 
one only needs to show that following equations have no 
solution with 9 = at far left and 9 = it at far right, 

^ = jy 2 9 HMsin9- Jsin0cos0(V0) 2 = 0. 

o9 o9 

^ = ^V-( S in^V0)-| = O. 

(2) 

Multiply the first equation by V0 and the second equa- 
tion by V(/>, then add up the two equations. One can 
show a tensor T satisfying V • T = with 



T =[-f + HMcos6+ ^{\V9\ 2 +s 
J(V0 <g> V6> + sin 2 9V(j) <g> V0), 



^|V^| 2 )]1- 



where 1 is 3 x 3 unit matrix. V0® V0 and V0<8> V</> are the 
usual dyadic products. The diagonal terms of T are just 
magentic Lagrangian density. If a HH DW exists with 
9 = in the far left and 9 — it in the far right, it requires 
-/(0, 0) + HM = -/(tt, (/)) - HAI that holds only for 
H = since f(0,4>) = f(ir,(f)). In other words, a static 
DW can only exist between two equal-energy-density do- 
mains. A HH DW in a nanowire under an external field 
must vary with time because two domains separated by 
the DW have different magnetic energy density. It should 
be clear that the above argument is only true for a HH 
DW in a homogeneous wire, but not valid with defect 
pinning that changes Eq. ((2|). Static DWs do exist in 
the presence of a weak field in reality because of pinning. 

What is the consequence of the non-existence of a 
static DW? A DW has to move when an external mag- 
netic field is applied along the nanowire as shown in 



Fig. 1. It is well known [12j| that a moving magnetiza- 
tion must dissipate its energy to its environments with a 
rate, ^ = eL ^- (dm / dt) 2 d 3 x , where m is the unit 

vector of M, a and 7 are the Gilbert damping constant 
and gyromagnetic ratio, respectively. Following the sim- 
ilar method in Ref. 13 for a Stoner particle, one can also 
show that the energy dissipation rate of a DW is related 
to the DW structure as 



dE cey 
~dt = ~ (1 + a 2 )M J_ c 



(M xH eff y d 3 x, (3) 



where H e f f — 



SF 
SM 



is the effective field. In regions I and 

II or inside a static DW (Fig. 1), M is parallel to H e ff, 
and no energy dissipation is possible there. The energy 
dissipation can only occur in the DW region when M is 
not parallel to H e f t . 

DW propagation and energy dissipation-For a magnetic 
nanowire in a static magnetic field, the dissipated energy 
must come from the magnetic energy released from the 
DW propagation. The total energy of the wire equals 
the sum of the energies of regions I, II, and III (Fig. 
1), E = Ej + En + Eu 1. Ei increases while En de- 
creases when the DW propagates to the right along the 
wire. The net energy change of region I plus II due to 
the DW propagation is d ^ gj j t gjI - > = — 2HMvA, where 
v is the DW propagating speed, and A is the area of 
wire cross section. This is the released Zeeman energy 
stored in the wire. The energy of region III should not 
change much because DW width A is finite, typically 
order of 10 ~ lOOnm. A DW cannot absorb or release 
too much energy, and can at most adjust temporarily en- 
ergy dissipation rate. In other words, dE d ' t 11 is either zero 
or fluctuates between positive and negative values with 
zero time-average. Since energy release from the mag- 
netic wire should be equal to the energy dissipated (to 
the environment), one has 
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2(1 + a 2 )HA J in 



(rh x i? e// ) 



d 3 x-\ 



dE 



in 



2HMA dt 
(4) 

Eq. ((4]) can serve as a proper definition of DW propaga- 
tion velocity that is completely defined by the instanta- 
neous DW structure. 

Velocity oscillation-A DW can have two possible types 
of motion under an external magnetic field. One is that a 
DW behaves like a rigid body propagating along the wire. 
This case occurs often at low field, and it is the basic 
assumption in Slonczewski model 11] and Walker's solu- 
tion for H < H~w ■ Obviously, both energy-dissipation 
and DW energy is time- independent, dE d ' t " = 0. Thus, 
the DW velocity should be time- independent. The other 
case is that a DW structure varies with time which oc- 
curs at large field H > Hw- In this case, DW structure 
deforms and DW precesses around the wire axis, expe- 
riencing different transverse magnetic anisotropy energy. 
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As a result, the DW width breathes periodically since it is 
defined by the balance between the magnetic anisotropy 
energy and the exchange energy. Thus, one should ex- 
pect both ^ff- and energy dissipation rate oscillate with 
time. According to Eq. ((H), DW velocity will oscillate. 
DW velocity may oscillate periodically or irregularly, de- 
pending on whether the ratio of precession period and 
breathe period is rational or irrational. Indeed, this os- 
cillation was observed in a recent experiment [8J] . How can 
one understand the wire- width dependence of DW veloc- 
ity? According to Eq. (j4|), the velocity is a functional of 
DW structure which is very sensitive to the wire width. 
For a very narrow wire, only transverse DW is possible 
while a vortex DW is preferred for a wide wire (large 
than DW width). Different vortexes yield different val- 
ues of |m x H e ff\, which in turn results in different DW 
propagation speed. 

Time averaged velocity is 



cry 



m x H. 



2(1 + <*)HA ./,„ V ") a ~ X > (5) 

where bar denotes time average. It says that the aver- 
aged velocity is proportional to the energy dissipation 
rate. The right-hand side of Eq. ([5]) is positive and non- 
zero since a time dependent DW requires m x H e ff =/= 0, 
implying a zero intrinsic critical field for DW propaga- 
tion. It is straightforward to show 
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, m x H. 



eff 



H sin 9 

M 



1 



M sin 9 



J dz 



sin 



dz 



df 



:d 2 



J d^'%- Js[n0cost 



(6) 



It is interesting to notice that two terms on the right 
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hand side of Eq. © are just [iff j ^ iM^e^e. 
In the case that there is no distortion on the DW plane, 



and [- 



8E1 



ie. d(f>/dz — and d 2 (j>/dz 2 — 0, then 
(m x H eff y = (H - H'q) 2 sin 2 9 - 



1 



df 



M sin 9 , 



(7) 
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FIG. 2: a) Symbols are from Ref. 8a. Solid line is the best fit 
to Eq. (8). Dash and dash-dot lines are fits to some modified 
expressions of Eq. (8). b) The time-averaged DW propa- 
gation speed versus the applied magnetic field for a biaxial 
magnetic nanowire. The cross are calculated velocities from 
Eq. J5), and the open circles are simulated average velocities. 
The solid curve is the fit to Eq. (8). Insets: the instantaneous 
DW speed calculated from Eq. © for H = 50Oe < H w (left) 
and H — lOOOOe > Hw (right), c) Symbols are numerical 
simulations of another set of parameters. Fhe dashed line is 
the linear fit for low field and solid lines are the fits to Eq. 
(8). d)-f) The dependence of a, Ho, and b on K' 2 defined in 
the text. Dash lines are linear fits, and solid line is quadratic 
fit. 



where H' = 



Alt 



</§j! — f^J ■ To reproduce the re- 
sult of Schryer and WalkeriQ for H < Hw, one sets 
|| = 0, H = 0, and J0 = || (similar to the origi- 
nal analysis of Schryer and Walker Q) so that H' = 0. 
Under the assumption, LLG Eq. gives a|| = 7_ff sin# 



andf = 



df 



-n—- — ~ -— • Thus , t 1 „ - Jr — — — sin 9, and 

M sin d(f> M sin aq) a ' 

the right-hand side of Eq (J7J) is just ^^-H 2 sin 2 9. Us- 
ing the definition of DW width A, / sin 2 9d 3 x = 2AA, 
then Eq. ([5]) reproduces the famous Walker's mobility 
coefficient (j, = 3— for H < Hw- For H > Hw, 4> varies 
periodically with time between and 27r, Substitute ex- 
pression in Eq. ((6]) into Eq. (|5|) and take the time av- 
erage, field dependent averaged velocity takes a form of 



v = aH — clq + a_i/ ' H that can be rewritten as 

v = a(H-H j 2 /H + b/H, (8) 

where a is proportional to the averaged DW width, H 
and b (and ao, a_i) depend on DW structure and mag- 
netic anisotropy. 

To demonstrate the goodness of Eq. ([5]) , we fit the ex- 
periment data (symbols in Fig. 2a) from Ref. 8a by the 
expression (solid line in Fig. 2a) with a — 2.98m/ (s-Oe), 
b = 5630e ■ m/s, and Hq = HOe. The experimental mo- 
bility at large field is measured to be 2.5m/ (s ■ Oe) that 
compares well with a — 2.98m/(s • Oe). According to 
our theory, b should be proportional to ajAK%/(l + a 2 ) 
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{K 2 , defined later, measures the transverse magnetic 
anisotropy (TMA)), where A is the time-averaged DW 
width. Using material parameters^ M — 860 x 10 3 A/m; 
J = 13 x lCT 12 J/m; K X M 2 = bOOJ/m 3 (defined later); 
a = 0.01, and measured low-field mobility /i = 25m/(s • 
Oe); and Walker breakdown field Hw = 40e, the DW 
width for H < Hw is about A = lAnm, and TMA con- 
stant to be K 2 M — 34 x 10 3 A/m from A = and 
Hw = aK 2 M, the ID result from Walker's original pa- 
per. It is known that DW width should vary as the DW 
precess around wire axis for H > Hw- Although exact 
value of the averaged width A is not known from the ex- 
periment, the theory requires its value to be about 16nm 
in order to obtain the fitting value of b = 5630e • m/s. 
I6nm is a fair value for DW width. The good agreement 
between Eq. JHJ) and experimental results is not acciden- 
tal. In fact, if one changes the fitting formula slightly to 
a(H-H ) 2 + b/H or a(H-H ) 2 + b, the best fits (dashed 
line and dashed-dot line in Fig. 2a) not only show poor 
agreement with experiment, but also give unreasonable 
fitting parameters. This proves that a good fitting is not 
due to three fitting parameters introduced. 

To further test the validity of Eq. (jHJ) and usefulness of 
both Eqs. ^ and ([5]) in evaluating the DW propagation 
speed from a DW structure, we carry out micromagnetic 
simulations on a strap wire of 4nm x 20nm x 3um whose 
magnetic energy density is F — —M ■ H + j4(VM) 2 — 
K x Ml + K 2 M 2 . We use OOMMF packageflj to find 
the DW structures and then use Eq. |5]| to obtain the 
average velocity or Eq. Q for instantaneous velocity. 
Fig. 2b is the simulation results for system parameters 
of Ki = 1Q- 7 N/A 2 , K 2 = 0.8 x IQ^N/A 2 , J = 4 x 
10- n J/m, M = 10 6 A/m, and a = 0.1. The calculated 
velocities are denoted by cross and numerical simulations 
are the open circles with their error bars smaller than 
the symbol sizes. The good overlap between the cross 
and open circles confirm the correctness of Eq. {5| . The 
V — H curve for H > Hw can be fit well by Eq. (jSJ. 
The insets are instantaneous DW propagation velocities 
for both H < H w and H > H w , by Eq. Q from the 
instantaneous DW structures obtained from OOMMF. It 
should be emphasized that both <f> and 6 are complicated 
functions of the time and the space. The left inset is the 
instantaneous DW speed at H = 50Oe < Hw, reaching 
its steady value in about Ins. The right inset is the 
instantaneous DW speed at H = lOOOOe > Hw, showing 
clearly an oscillation. They confirm that the theory is 
capable of capturing all the features of DW propagation. 

Figure 2c is another set of OOMMF simulations on 
a wire of Anm x 16nm x 4/im with system parame- 
ters of M = 500k A/m for the saturation magnetization, 
J = 20 x 10~ 12 J/ra, the axial crystalline anisotropy con- 
stant Ki = 8 x 10~ 7 N/A 2 , and a = 0.1. The symbols are 
simulation results, and solid lines are the fits to Eq. ©. 
The perfect agreement demonstrate the correctness of 
Eq. ([5|). Fig. 2d-2f show how a, b and Hq depend on the 



normalized TMA constant K' 2 = K 2 + 3.4 x 10~ 7 N/A 2 . 
This normalized parameter comes from the demagnetiza- 
tion effect that generates an extra magnetic anisotropic 
energy (D 1 -D 2 )M^ + (D 3 -D 2 )M^+D 2 , where A is the 
demagnetization factor along i-th axis. Using the theory 
in Ref. 15, one can find D 1 - D 2 to be 3.4 x 10" 7 A/A 2 
for our geometry, a and Hq are linear in K 2 while b is 
quadratic in K 2 , in agreement with Eq. (7). It should 
be pointed out that velocity expressions in both Refs. 
9 and 10 cannot fit either the experimental curve (Fig. 
2a) or simulations (Figs. 2b and 2c). The correctness of 
result Eq. (fj)) depends only on the LLG equation and 
the general energy expression of Eq. ([1]) . It does not de- 
pend on the details of a DW structure whether they are 
transverse or vortex like. In this sense, our result is very 
general and robust, and it is applicable to an arbitrary 
magnetic wire. 

In conclusion, a proper definition of DW propagation 
velocity is obtained, and a velocity-field formula for high 
field (above the Walker breakdown field) is proposed. 
This new formula agrees well with both experiments and 
numerical simulations. Furthermore, a global picture of 
DW propagation in a nanowire driven by a magnetic field 
is revealed: A static DW cannot exist in a homogeneous 
magnetic nanowire when an external magnetic field is ap- 
plied. Thus, a DW must vary with time under a static 
magnetic field. A moving DW must dissipate energy due 
to the Gilbert damping. As a result, the wire has to 
release its Zeeman energy through the DW propagation 
along the field direction. The DW propagation speed is 
proportional to the energy dissipation rate that is deter- 
mined by the DW structure. An oscillatory DW motion, 
either the precession around the wire axis or the breath 
of DW width, should lead to the speed oscillation. The 
observed negative differential mobility is due to the tran- 
sition of a DW from a high energy dissipation state to a 
low energy dissipation one. 
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